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A theory is developed for the T — Mott-Hubbard insulating phases of the d°° Hubbard model 
at j-filling: including both the antiferromagnetic (AF) and paramagnetic (P) insulators. Local 
moments are introduced explicitly from the outset, enabling ready identification of the dominant 
low energy scales for insulating spin-flip excitations. Dynamical coupling of single-particle processes 
to the spin-flip excitations leads to a renormalized self-consistent description of the single-particle 
propagators that is shown to be asymptotically exact in strong coupling, for both the AF and P 
phases. For the AF case, the resultant theory is applicable over the entire [/-range, and is discussed 
in some detail. For the P phase, we consider in particular the destruction of the Mott insulator, 
the resultant critical behaviour of which is found to stem inherently from proper inclusion of the 
spin-flip excitations. 
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1. INTRODUCTION 

Since its inception more than thirty years ago pi, the 
Hubhard model has become the canonical model of in- 
teracting fermions on a lattice. Although possibly the 
simplest model to describe competition between electron 
itinerancy and localization, with attendant implications 
for a host of physical phenomena from magnetism to 
metal-insulator transitions, its simplicity is superficial 
and an exact solution exists only for d = 1 dimension 

Recently, Metzner and Vollhardt |3j have pointed to 
the importance of the opposite extreme, d = oo. In 
suppressing spatial fluctuations, the many-body prob- 
lem here simplifies considerably, reducing to a dynamical 
single-site mean-field problem. Motivated in part by the 
expectation that an understanding of the d = oo limit 
will serve as a starting point for systematic investigation 
of finite dimensions, and by the knowledge that some 
important vestiges of finite-d behaviour remain inherent 
in the d — oo limit, intense study of the ^-filled d = oo 
Hubbard model on bipartite lattices has since ensued; for 
recent detailed reviews, see Refs f|,[|,|]]. 

The true ground state of the model is an antiferromag- 
net (AF) for all interaction strengths U > 0. One aim of 
the present work Q is to develop a theory for the d = oo 
AF which, in contrast to previous theories for the AF 
phase pPJlC}], is reliable over the entire £/-range, and in 
particular becomes exact in the U — > oo strong coupling 
limit both at i-filling where the Hubbard model maps 



onto the AF Heisenberg model, and in the one-hole see- 




on the d = oo 



tor where it reduces to the t-J model 

The majority of previous work [| 
Hubbard model has focused on the paramagnetic (P) 
phase that results, even for T = 0, simply by neglecting 
the magnetic ordering (or suppressing it via frustration 
H ) . One highlight of this work has been the emergence 
of a detailed description of the Mott metal-insulator tran- 
sition, although here too the picture is not complete: 
for example, a firm understanding of the mechanism by 



which the T = Mott insulating solution is destroyed, 
and even whether it is continuous or first-order, remains 
elusive f^,[f2|]. A second aim of this paper is to focus on 
the insulating state of the P phase, and to develop a the- 
ory for it on a footing essentially identical to that for the 
AF, which likewise becomes exact in strong coupling and 
which permits an analysis of the destruction of the Mott 
insulator. 

In seeking to develop a 'unified' description of the AF 
and P insulating phases, we adopt a rather different ap- 
proach to that taken in previous work by intro- 
ducing explicitly, and from the outset, the notion of site 
local moments. To this end we consider first a conven- 
tional T = mean-field approach to the problem in the 
form of unrestricted Hartree-Fock (UHF), together with a 
random phase approximation (RPA) for transverse spin 
excitations of the mean-field state. Despite the limita- 
tions of such an approach per se, its importance resides 
in enabling identification of key low energy scales for in- 
sulating spin-flip excitations. Since spatial fluctuations 
are suppressed for d — oo, the low energy spin- flip ex- 
citations are found to be Ising-like and (for each phase) 
characterized by a single scale, u/ s . This has a simple 
physical interpretation. For the AF, lu s — u p (U) is es- 
sentially just the energy cost of flipping a spin in the Neel 
ordered background; the ubiquity of antiferromagnetism 
for all U > leads naturally to u> p > for all U, with 
ujp ~ 1/U as U — > oo as one expects in the Heisenberg 
limit. For the P phase by contrast, where magnetic or- 
dering is absent, the fact that a given spin is equally likely 
to be surrounded by f- or J,-spins and thus (for d = oo) 
has as many j- as |-spin neighbours, ensures that the 
corresponding spin-flip energy cost u> s = for all U in 
the insulating state. 

Identification of the low energy spin-flip scales, while 
crucial to the present work, is preliminary: to transcend 
the limitations of the conventional 'static' mean-field ap- 
proach, single-particle processes must subsequently be 
coupled dynamically to the transverse spin-flip excita- 
tions. It is this which, in leading as we shall describe to 
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a self-consistent description of the single-particle Green 
functions, enables the aims outlined in the preceding 
paragraphs to be achieved. 

The Hubbard Hamiltonian, in standard notation, is 

H = -t } C^Cja +U22 n it n il -t = U/ \[2Z 
(ij)a i 

(1.1) 

with the {if) sum over nearest neighbour sites on a bi- 
partite lattice of coordination number Z: a Bethe lat- 
tice (on which in practice we shall largely focus), or a 
d-dimcnsional hypercube. To ensure a non-trivial limit 
as d — > oo ||, the hopping is scaled as t = t*/y/2Z. The 
paper is organised as follows. UHF+RPA, and the spin- 
flip scales referred to above, are discussed in §2. Em- 
phasis is also given here to simple physical arguments 
which, in highlighting both the deficiencies and virtues 
of UHF+RPA, indicate what is required to go beyond it; 
particular attention being given in this regard to UHF for 
the P phase, in view of its close relation to the early work 
of Hubbard ftf] and the Falicov-Kimball model |j. 

Dynamical coupling of single-particle processes to the 
transverse spin excitations is considered in §3, leading 
(§3.1) to a renormalized self-consistent approximation for 
the (T = 0) single-particle Green functions upon which 
we subsequently concentrate. In §3.2 the strong cou- 
pling behaviour is examined analytically, and shown to 
be asymptotically exact for both the P and AF phases. 
Results are given in §3.3, focusing in particular on single- 
particle spectra for the AF from strong to weak coupling, 
and on a discussion of the localization characteristics of 
the single-particle excitations -the latter being quite sub- 
tle, and pointing to the delicacy of the limit U — > oo for 
the AF phase. For the P insulator, single-particle spec- 
tra are discussed briefly in §3.3, before considering the 
destruction of the Mott insulating solution in §4. The 
single-particle gap is found to close continuously, with 
an exponent v = 1, at a critical U c — 3.41t». The ori- 
gins of this behaviour are found to stem from inclusion of 
the io s = spin-flip scale in the interaction self-energies, 
pointing to the importance of such throughout the en- 
tire insulating regime, and not solely in obtaining the 
exact strong coupling limit. The results of §4 are in good 
agreement with recent numerical work p2| , as discussed 
in §5. 



2. CONVENTIONAL MEAN-FIELD APPROACH 

We focus on the zero temperature single-particle Green 
functions, defined by 

Gu, a = -ma^cl}) = G± a {t) + Gr.^t) (2.1) 

(for the site diagonal element); and separated for later 
purposes into retarded (+, t > 0) and advanced (— ,t < 0) 



components. The essential feature of d°° is that the cor- 
responding interaction self-energy is site-diagonal |§,[l5) , 



(0) 



here, and throughout, Co denotes 



frequency relative to the Fermi level, viz Co — lo — U/2. 
Gn- a {uj) may be written as 



G 



i;o-(w) = [u> - Si CT (w) - S ia (uj)] 1 (2.2a) 

where Si„ is the 'medium' self-energy — which alone sur- 
vives in the non-interacting limit — expressing hopping of 
CT-spin electrons to neighbouring sites. Simple application 
of Feenberg's renormalized perturbation theory (jHj 17 
shows that, for d = oo but regardless of lattice type, Si a 
is a functional solely of the {Gjj.a}. The functional de- 
pendence is particularly simple for the Bethe lattice (BL) 
on which we concentrate, namely 
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{Co) 



(2.2b) 



with Uj = t se /\'2Z the nearest neighbour hopping ele- 
ment. Note that this is quite general; no assumption has 
been made about magnetic ordering or otherwise. 

We consider now a conventional mean-field approach 
to the single-particle Green functions. 



2.1 UHF 

For both the AF and P phases, a Hartree-Fock ap- 
proximation — by which we emphasize is here meant spin 
unrestricted Hartree-Fock (UHF) — is the simplest non- 
trivial mean field approximation, in which the notion of 
site local moments (/ij), regarded as the first effect of 
electron interactions, enters from the outset. In the AF 
case, the local moments are naturally ordered in an A/B 
2-sublattice Neel state, with Uj — ±|/x| for site i in the 
A/B sublattice respectively [|18|. For the P phase by con- 
trast, the local moments are randomly oriented: a site is 
equally likely to be A- type as B-type [jl!|. In either case 
the essential — and limiting — feature of UHF is that it 
is a static approximation, with solely elastic scattering 
of electrons and w-independent interaction self-energies 
approximated by 



y0 



(2.3) 



For the AF phase, the UHF Green functions (G° . a 



G° aa with a =A or B) follow from Eqs {23 



BL as 



G Aa {Co) =[£ + §[/| M | 
Gl„{u>) =[u>-fC/| M | 



¥IgL[Co)\ 
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for the 

: AF 
(2.4) 



where the 'medium' self-energy part reflects the 2- 
sublattice structure of the Neel state. Eqs (2.4) are a 
closed set, with the UHF local moment |^| = found 
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self-consistently via the usual gap equation (see e.g. p0[), 
which may be written formally as 



iMol = 



du (D A f(u)) - D a1 (lu)) 



(2.5) 



in terms of the corresponding spectral densities. And the 
total Green function is given by 



(2.6) 



such that D°(uj) = — 7r~ 1 sgn((I))ImG°(£D) gives the total 
single-particle spectrum. 

For the P phase by contrast, 



Gl a {u>) =[Q + !LU\n\-UlG° 
G% a {Q) =[cD-fC/|M|-|^G (a))]- 1 



P (2.7) 



The sole difference to E q. ( p.4[ ) occurs in the medium 
self-energy (see Eq. ( |2.2| b)), since the nearest neighbours 
to any site are equally likely to be A- or B-type sites. 
Eqs are a closed set for G°(uj) and the G° a(T (Q) in 

the P ph ase; the UHF local moment is again found from 
Eq. (£f). 

For either phase there are two basic symmetries, viz 



D Aa {{b) = Di_„{Q) 



(2.8a) 
(2.8b) 



reflecting the | / J.-spin symmetry (G ArT (uj) = G®_ a (oj)) 
and particle-hole symmetry (G Aa (uj) = — G° } _„(— uj)) re- 
spectively; and note therefore from Eq. ( |2.6| ) that G°(Q) 
is naturally independent of the spin, a. 

We add further that UHF yields the correct atomic 
limit (where kip | = 1) for either phase, as is clear from 
Eqs (|j,|jU2j) with i* = 0. 



A. Antiferromagnet 

UHF for the AF has been widely studied since the 
early work of Penn jTsj. Here we mention only that for 
any d > 1 the exact ground state of the i-filled Hub- 
bard model on a bipartite lattice is an AF insulator for 
all U > 0, and this is qualitatively well captured at UHF 
level: for all U > 0, the mean-field ground state is a 
2-sublattice Neel AF, with a gap in the single-particle 
spectrum D°(uj) given by A(U) = U\fi\; Fig. 1 shows 
D°(u>) at U/U = 4 for the d°° BL. 

The deficiencies of UHF are however most clearly seen 
in strong coupling, U — > oo, where for the AF the single- 
particle spectrum reduces to D°(uj) = ^[S(u>+^) + 6(u> — 
¥)] — as for the atomic limit, t„ = 0. The physical origin 
of this is simple: consider for example the upper Hubbard 
band in strong coupling, and imagine adding an |-spin 
electron to a site (B-type) already occupied by a .[-spin. 
Since UHF is an independent (albeit interacting) electron 
approximation, only the added [-spin can potentially hop 



to nearest neighbour (NN) sites. But it cannot do so in 
the strong coupling limit, since for the AF all NN's to the 
!-spin B-site are |-spins (A- type). The added T-spin thus 
effectively 'sees' the |-spin site as an isolated site, hence 
the emergence of atomic limit behaviour as the strong 
coupling limit at UHF level. But while physically trans- 
parent, this behaviour is wrong. In strong coupling, and 
for the 1-hole sector appropriate to the lower Hubbard 
band (or 1-doublon sector for the upper Hubbard band), 
the Hubbard model maps onto the t-J model Jll]] 



H tJ = -t 



E4 



\Joo ^ ' (&i ' Sj 



irnnj) (2.9) 



where the hole moves in a restricted subspace of no dou- 
bly occupied sites (cj^ = ct ff (l — ni_ a )); and in the fluctu- 
ating spin background provided by the Heisenberg part of 
HtJ, with NN exchange coupling Joo = At 2 /U. Although 
it is exact in the atomic limit, UHF by itself can evidently 
say essentially nothing about the strong coupling limit. 



B. Paramagnet 

UHF for the T — paramagnetic phase warrants sep- 
arate discussion, in part because of its very close rela- 
tion to two other well known approaches. The first is 
that due to Hubbard with 'spin-disorder scattering' 
only. This is often called the Hubbard III (HIII) approx- 
imation, and we here refer to it thus (noting that 'reso- 
nant broadening' contributions are additionally included 
in Ref. (l3)). HIII is equivalent to UHF, but with a satu- 
rated local moment. Thus, with = 1, the resultant cu- 
bic equa tion for G°(oj) on the d°° BL, obtained from Eqs 
( [Hi 2.7) above, coincides precisely with the HIII approx- 
imation for any U ; see eg Eq. (34) of Ref. Although 
Hubbard's original formalism is very different, its physi- 
cal content is that of a static approximation to an alloy 
analogy description |p2fl ; a close relationship to UHF is 
thus to be expected. 

The second connection is to the Falicov-Kimball (FK) 
model |i~4|], a simplified version of the Hubbard model 
in which electrons of only one spin type are mobile, and 
which for d°° is exactly soluble [ p3p^ , |25] | . For the para- 
magnetic phase of the FK model the single particle Green 
function reduces precisely to that of HIII for any U [^2| , 
ie to the above-mentioned cubic for G°(uj) on the c?°° BL; 
see eg Eqs (7.3) and (4.4) of Ref. ||. 

For the P phase, Fig. 1 shows the UHF D°(w) for 
U/t* = 4 on the d°° BL, contrasted to its AF counter- 
part. Ordering or otherwise of the preformed local mo- 
ments clearly has a significant effect on the spectra. In 
the AF ordered case, for example, the interior edges of 
the Hubbard bands have characteristic square-root diver- 
gences, while for the P phase all band edges vanish with 
square-root behaviour. More significantly, while the 2- 
sublattice structure of the Neel ordered state ensures a 
band gap A = U\fi\ for all U > 0, the single-particle UHF 
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gap vanishes in the P phase at a critical U c ~ 1.9f* given 
by U c \[i(U c )\ = \/2i* — or correspondingly U c — \/2t r 
for HIII/FK — signalling an insulator-metal transition. 
UHF/HIII fails of course in the metallic phase, there be- 
ing no well-defined Fermi surface or quasiparticles |^2|,^6| . 
This is inevitable for any inherently static approximation 
with a frequency-independent self-energy E Q(7 , since the 
essence of Fermi liquid behaviour is the inelasticity of 
electron scattering near the Fermi level, Q — (^?J . 

However, even in the P insulating phase of interest 
here, UHF/HIII is deficient. As for the AF this is seen 
most clearly in strong coupling, U — > oo, where although 
the centres of the Hubbard bands are separated by U, 
each has a non- vanishing width. In the strong coupling P 
phase, and for the one hole (doublon) sector correspond- 
ing to the lower (upper) Hubbard band, the Hubbard 
model maps onto the t-J model (Eq. (2.E)) in a random 
spin background; and the exact full bandwidth of either 
band is given for the d°° BL by [t28||2S 



U 



(2.10a) 



Note that this is also the single-particle bandwidth in 
the other extreme of the non-interacting limit, reflect- 
ing physically that in strong coupling the hole/doublon 
behaves essentially as a free particle Q. 

In contrast, the strong coupling bandwidth at 
UHF/HIII level is 



2U 



UHF/HIII 



(2.10b) 



UHF or HIII does not therefore give the exact strong 
coupling limit for the Hubbard model, contrary to what 
has been suggested recently but, as is clear from the 
above discussion, gives instead the strong coupling limit 
of the FK model |£5). The physical origin of Eq. fl2.10|b) 
is however both simple and revealing. Consider again 
the upper Hubbard band in strong coupling, and imagine 
adding an f-spin electron to a B-type f-spin site. Within 
a static approximation such as UHF/HIII, only the added 
f-spin can hop; and it can do so in the first instance to 
any of J. -spin NN's (B-type sites) only — the effective co- 
ordination number for which is Z c g = ^Z. Since Z e g 
for the propagating "f-spin electron is reduced by a factor 
of 2 below the full coordination number, and since the 
bandwidth of the d°° BL is proportional to V-^off > the 
strong coupling UHF/HIII width is thus diminished by 
y/2 from t he co rresponding non-interacting value which, 
as in Eq. (|2.10|a ) above, is also the exact strong coupling 
limit; Eq. (|2.10|b) thus results. 



The distinction between Eqs (2.10a) and (2.10b) is 
however qualitative, and not solely a matter of de- 
gree, reflecting the need to take seriously — even in 
strong coupling — the correlated dynamics of the elec- 
trons. Whenever, say, an f-spin electron is added to 
a site occupied by a f-spin electron, the added f-spin 
can indeed propagate in the P phase, scattering elasti- 
cally off successive neighbouring f-spins; and as sketched 



above this is well captured at UHF/HIII level. But, hav- 
ing added the f-spin to a f-spin site, the latter can itself 
clearly hop off the site — to a neighbouring f-spin site — 
leaving behind it a spin-flip on the original site. The en- 
ergy cost for the spin-flip is zero, since we are considering 
the P phase where, for d°°, a given spin is equally likely to 
be surrounded by f or f spins and has as many f as f-spin 
neighbours (whence there is no 'exchange penalty' for a 
spin- flip). Thus, whether the added f-spin or the f-spin 
already present hops off the site, the initially create d dou- 
blon propagates as a free particle p0[ ]; Eq. ( 2.10| a) thus 
results. To describe correctly the electron dynamics, both 
types of process above — and therefore the interference 
between them — must be included. A static approxima- 
tion such as UHF/HIII cannot handle this, since such 
dynamics reside in the frequency dependence of the full 
interaction self-energy E^u)), as considered in §3. 



2.2 RPA 



In contrast to single-particle spectra — probing states 
one hole or particle away from ^-filling — RPA probes 
fluctuations about the mean-field state, and thus ex- 
citations of the ^-filled state itself. For the insulat- 
ing phases, with a gap to charge excitations, transverse 
spin excitations are of lowest energy. These are re- 
flected in the transverse spin polarization propagators 



TJ±-(t) = i(T{S+(t)S~}) and IE+ft), given within RPA 
by 



°n+- (uj)[i - u °n + -( w )]^ 1 (2.iia) 



Sij and u nJ is 



where [n+-(w)]« = U+~{u), [% _ 
t he p ure UHF transverse spin polarization bubble. Eq 
fl2.1l[ a) leads directly to a familiar diagrammatic 'bub- 
ble sum'. Alternatively, since the interaction is solely 
on-site, this may be recast as a 'ladder sum' of repeated 
particle-hole interactions in the transverse spin channel, 
as shown in Fig. 2; bare UHF propagators are denoted by 
solid lines, and the on-site interactions (conserving spin 
at each vertex end) by wiggly lines. From the basic sym- 
metries (Eq. (p|), it follows that U^ + (lu) = H±-(-w) 
for i = A or B; and U+~ {lj) = U~+(uj). 

For finite-<i, intermediate sites in the ladder sum for 
(Fig. 2) are in general different from i. But since 
°Hij ~ 0(d~ rn ) for sites i and j m th nearest neighbours, 
all intermediate sites in are equal to i for d°°; ie 
i = i\ = %2 = ■ ■ ■ = i whence Hj _ (or n^ + ) is purely 
algebraic, viz 



n+-H= °n+-( w )/[i-c/°n+-H] 



(2.11b) 



The spectral density of transverse spin excitations is 
naturally reflected in Imllj~(a;), as now considered for 
the AF and P phases. 
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A. AF phase 

For the AF, Fig. 3a shows ImU^w) at U/U = 4 
for the d°° Bethe lattice. Two distinct features are ap- 
parent: a low frequency spin- flip pole (discussed below), 
and a high energy Stoner-like band. The latter consists 
simply of weakly renormalized Hartree-Fock excitations 
across the gap in the mean-field single-particle spectrum. 
Spectral density for the Stoner bands does not therefore 
begin until precisely \u>\ = U\fi\ (see Fig. 1), and their 
maximum density occurs for \u>\ ~ U. This is as found 
also for fmite-ei, see eg Ref [p2| . 

The central fe ature in Fig. 3a is a law-u> pole at lo p , 
located via Eq. ( 2.11 b) from U ~ (w p ) = 1, and oc- 
curring for all U > (Fig. 3a, inset). This is the sole 
remnant, for (i 00 , of the spin wave-like component of the 
transverse spin spectrum studied recently at RPA level 
p2|; and which, for finite d and any U > 0, is natu- 
rally gapless. Physically, the single spin- flip pole at lj p 
reflects the general suppression of spatial fluctuations for 
d°° : it corresponds simply to the energy cost of flipping a 
spin in the AF background. This is particularly clear in 
strong coupling, where the Stoner bands are eliminated 
entirely. Here, as is well known 20 33], the RPA trans- 
verse spin spectrum reduces (for any a) to the linear spin 
wave spectrum of the nearest neighbour AF Heisenberg 
model, with exchange coupling = 4t 2 /U = 2t^/ZU, 
onto which the ^-filled Hubbard model maps rigorously. 
And for d°° it is straightforward to show that the resul- 
tant linear spin wave spectrum collapses to an Ising-likc 
spin-flip pole at = ZJoc/2 = t\jU . Further, since 
linear spin wave theory for the Heisenberg model is exact 
for d°° Q, it follows that in strong coupling UHF+RPA 
gives the exact spin excitation spectrum of the ^-filled 
Hubbard model. 

The occurrence of the single w p -pole is robust to fur- 
ther renormalization of particle- hole lines in Y[^~(ui), as 
discussed in §3.3. We stress further that to capture it 
requires the full ladder sum of repeated p-h interactions 
shown in Fig. 2: retention solely of the 'bare' polariza- 
tion bubble diagram will clearly not suffice. 

The necessity of including the AF spin-flip scale will 
be evident when discussing the T — single-particle 
spectra, §3.2,3. Here, we illustrate briefly its impor- 
tance at finite temperature, as reflected in the Neel tem- 
perature Tm{U). Molecular field theory is exact for 
the Heisenberg model in d°° pa] ; thus, in strong cou- 
pling, Tjv = ZJoo/A — At finite U, we expect 
Tjv(J7) ~ \tu p {U) to yield a good estimate of the Neel 
temperature in a U -regime where thermal properties are 
dominated by the low- lying spin- flip excitations. Jarrell 
and Pruschke |36|j3^ ] have obtained the finite-T phase di- 
agram for the d°° hypercubic lattice via quantum Monte 
Carlo (QMC). The thermal paramagnetic phase above 
Tn(U) is found to be metallic for /7/i* 3 and insulating 
for U/t* > 3 (with a small 'crossover' regime); it is thus 



borne out. Fig. 4 shows the QMC T N (U) for the d°° hy- 
percubic lattice, together with the corresponding ^u) p (U) 
and the exact strong coupling asymptote Tat = ^w^ . 
The QMC Neel temperature is indeed well described by 
t}U) p (U) down to U/t* ~ 3. 

B. P phase 



For the T = P phase, Fig. 3b shows Imn+~ (oj) at 
U /i* = 4 for the d°° Bethe lattice. Compared to its 
AF counterpart (Fig. 3a) the key difference is that the 
spin- flip pole occurs at oj = 0, reflecting the fact that 
the energy cost for a spin-flip is zero in the paramagnetic 
insulator, as argued physically in §2. lb. The formal ori- 
gin of this at RPA level is seen readily by noting that the 
bare transverse spin polarization bubble (Fig. 2, diagram 
(a)) is given by 



— G Ai (Q')G%(u;'-u;) . (2.12) 



From Eq. (|2.4|) f° r the UHF Green functions, using 



G° Al = G° bV it follows that 



G° Al (u,)G%(u,) 



u\/A 



{G%(u>)-G Al (u)) . (2.13) 



in the latter region that we expect Tj 



N 



This is 



Hence, using the spectral representation of G A(T (ui), 

°H+-( W = 0) = -|- f du [D° Ar (Q) - D° Al (u>)] . 

(2.14) 

Sinc e the UHF local moment = \(1q\ is given by E q. 
(U), n+-(w = 0) = l/U; and thus from Eq. (pTl| b) 
the RPA (oj) has a spin-flip pole at ui = Q. 

Note again, as for the AF case, that the full ladder sum 
of particle-hole interactions in the transverse spin channel 
is required to capture the u) — spin-flip pole. Further, 
although we have shown explicitly its existence within 
RPA, the occurrence of the zero-frequency spin-flip scale 
is naturally a general feature of the d°° paramagnetic in- 
sulating phase where, locally, the ground state is a dou- 
bly degenerate local moment (as for the single-impurity 
Anderson model embedded in an insulating host) ||. 

For both phases, the evident virtues of the RPA for 
excitations of the i-filled state contrast sharply with the 
deficiencies of the single-particle spectra at UHF level, 
§2.1. This itself hints at what is necessary to describe the 
single-particle spectra successfully: single-particle pro- 
cesses must be coupled dynamically to the transverse spin 
excitations, reflected in the frequency dependence of the 
self-energy. This is now considered. 
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3. GREEN FUNCTIONS 

It is helpful to separate the full interaction self-energies 
S acr (<D) as 



£ ACT (£>) 

S B( r(&) 



--L%| + £ a<t i 



0) 



-t/|/z| + £ BtT (w) 



(3.1) 



where £ atT (w) (a =A or B) excludes the first-order UHF- 
type contribution, and contai ns the d ynamics on which 
we want to focus. From Eqs (2/2, 23) for the Bethe lat- 
tice, the exact site-diagonal Green functions are thus 
given formally by 

G A(T (cD) = [Q+ a -V\ii\ - S Aa (u) - E^Cw)]- 1 (3.2a) 

G Ba (u) = [u> - ^U\n\ - S Ba {u>) - Huaiu))]- 1 . (3.2b) 

Here, the medium self-energy is given for the AF and P 
phases by 



S a A&) = { 1.2 



hlG &a (u>) ■ AF 
2 t 2 Mu>) : P 



(3.2c) 



where the site index a =B or A for a —A or B respec- 
tively; and 



G(Q) 



1 



[G Aa {u>) + G Ba (u) 



(3.3) 



is the total Green function. 

As at UHF level, f / |-spin symmetry and particle- 
hole symmetry for the corresponding spectral densities 
imply 



D A<y {Cb) = £> B _ CT (u>) 
= D A - a (-w) 



(3.4a) 
(3.4b) 



respectively. For the associated Green functions 
G aa (Cj) — G~£ a (ui)) + G~ a (w), a Hilbert transform of Eqs. 
( |3.4|) gives directly 



Thus, from Eq. (|J), 

G(u>) = -G(-w) ; 



while from Eq. (3.2,3.5) 



S B(T (w) = E A _ CT (a}) 
= -£ ACT (-£) 



(3.5a) 
(3.5b) 

(3.5c) 



(3.6a) 
(3.6b) 



and likewise for the S Q(T 's. Eqs ( |3.5| a) with ( |3.3| ) shows 
also that G(uj) is correctly independent of spin. 

The symmetries reflected in Eqs ( 3 . 5 , |3 . 6| ) play an im- 
portant role in the following analysis. For the P phase, 



note also the physical interpretation of Eq. (3.3) for 
G(d>): viewing the paramagnet in terms of randomly ori- 
ented local moments, where a site is equally likel y to 
likely to be A-type as B-type, we can consider Eq. (3.3) 
as a configurationally averaged Green function. This is a 
natural alloy analogy interpretation but, unlike the static 
approximation to such inherent in UHF or HIII, it is for- 
mally exact since no approximation to the interaction 
self-energies has thus far been made. 



3.1 Self-consistent renormalization 

Our aim now is to develop a specific approximation 
to the self-energy which in particular (a) becomes exact 
in strong coupling, ensuring thereby a controlled limit; 
and (b) is constructed in renormalized form, enabling a 
self-consistent solution for the single-particle Green func- 
tions. 

A relevant diagram contributing to is shown in 
Fig. 5, employing the same diagrammatic notation as 
Fig. 2. Using deliberately a strong coupling terminol- 
ogy, its physical interpretation is as follows (with t > 
for convenience): at t = a (a =)f-spin electron, say, 
is added to site i, thus creating a 'doublon'; at ii > 
the (— cr=)|-spin electron already present on site i hops 
from i to j, and at t% > t\ an "["-spin hops from j to k; 
the entire path is then retraced. The diagram thus de- 
scribes motion of the doublon (or hole for t < 0) from 
i j — > k via a correlated sequence of alternating spin 
hops, creating behind it a string of flipped spins. All 
ladder interactions of the resultant on-site particle-hole 
pair — which reflect the on-site spin-flip created by mo- 
tion of the doublon/hole — are shown explicitly for site i 
in Fig. 5; from which it is seen that their sum is exactly 
C/ 2 n^+, with Ur+(Lu) (= n+ - (-w)) the RPA transverse 
spin propagator discussed in §2.2 (cf Fig. 2). 

It is precisely correlated dynamics of the sort exempli- 
fied by Fig. 5 that we seek to include and generalize in 
the frequency-dependent E Q(T (a)). To this end we first de- 
fine an undressed (or self-consistent host) Green function 
by 



(3.7) 



This is shown dia gram matically in Fig. 6(a), as obtained 
simply from Eq. (3.7) using the Dyson equation for the 
full Green function Gu y(7 expressed in terms of the UHF 
propagators and self-energy insertions. As seen from the 
figure, the implicit sum over intermediate sites j, k etc. 
is thus restricted to exclude site i itself (unlike the full 
Gu ]CT where the site sums are free). While including all 
interactions on sites j ^ i, Gu-a- thus excludes all interac- 
tions on site i beyond the simple first-order UHF contri- 
bution to Sj CT . The latter is of course subsumed into the 
UHF Green functions (as in §2.1), which here constitute 
the 'bare' propagator s; an d in this important sense the 
above definition, Eq. (3.7), of the host propagator differs 
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from that of eg Refs |y,g§g7|3|j3S|j40[ (which would be 
recovered if we set the site local moment \fi\ = 0). 

To generalize the processes contained in Fig. 5, we 
renormalized the self-energy as shown in Fig. 6(b), re- 
placing the — CT-spin particle lines connecting the starred 
vertices i in Fig. 5 by the self-consistent host Green func- 
tion Qu-.-a] the infinite set of diagrams thus retained in 
Yii a follows simply by direct iteration of Fig. 6(b) us- 
ing Fig. 6(a) for Qa - a . This renormalization is adopted 
for the following reasons, (i) It ensures that an on-site 
spin- flip occurs only when the doublon/hole hops off a 
site, and that its outward path is self-avoiding. With ref- 
erence to Fig. 5 for example, site j ^ i is guaranteed, 
likewise k ^ j; while terms with k = i vanish for 
being at least 0(l/d) since G° ;(T ~ 0(d~ m / 2 ) for sites i 
and j m th nearest neighbours, (ii) In addition, the resul- 
tant site restrictions further prevent the need to include 
a class of partially cancelling exchange diagrams, as il- 
lustrated simply in Fig. 7 (where a sum over j i is 
implicit). Since j ^ i is guaranteed, the exchange dia- 
gram (Fig. 7(b)) is at least 0(1/ d) and thus vanishes 
for d°°, while the 'direct' diagram (Fig. 7(a)) is 0(1). 
If, however, j = i was included in the direct diagram, its 
exchange counterpart would also be O(l) and would thus 
need to be retained. 

Our basic approximation to Tn a (Cj) is thus Fig. 6(b), 
namely 



S AT (w) = U 2 



(3.8a) 



the remaining S Qff 's follow by symmetry, Eq. ( p.6[ ). 
From §2.2 the RPA n"+(0) (= II+r(-f2)) may be sep- 
arated into the spin-flip pole contribution, Q[Q + lo s — 
iri]" 1 (with pole- weight Q), plus the Stoner contribution; 
may be cast as 



whence Eq. (3 

£ AT (£)=Q[/^(a). 
with spin-flip frequency: 

10 



v^Stoner y,~ 



AF 
P 



(Q) (3.8b) 



(3.8c) 



By symme try, <? A j.(f2>) = Gbi(&) using Eqs (3.5-7); and 
from Eqs ([0> 2 ) : 



£a|(£>) 



iJ7|/i| 
|C/| M | 



¥lG{Co)}- x 



: AF 
: P 



(3.9) 



From Eq. (3.8a), £ A f(d;) is thus a functional of the 
Green functions, the basic equations for which (Eqs ([3.2])) 
must therefore be solved self-consistently, as now de- 
scribed. 



3.2 Strong coupling 

We consider first the behaviour in strong coupling, 
as this can be extr acted analytically. Since = 1 — 
0(tl/U 2 ), Eqs (|3.2|) reduce in strong coupling to 



G aT (w) = [u - 5 aT (w) - EatM] 1 
G BT H = [u-U- S b] -(lu) - Sbt(w)]" 1 



(3.10a) 
(3.10b) 



The interaction self-energy, Eq. (3.5), likewise simplifies 
in strong coupling, since the Stoner contribution vanishes 
(see^ §2.2) and the pole-weight Q — > 1. Hence from Eq. 

C/~ 2 £ aT (cj) = g^(uj + u> s ) = g~j(u> + u s ) . (3.11) 

This may be reduced further, noting that Q~Auj) is given 
by 

GM = [ U/2 *». PbtM . (3.12) 

as a one-sided Hilbert transform of the correspond- 
ing lower Hubbard band spectral density, P B t(^i) = 
— 7r _1 sgn(a;i — U/2)lmQ B ^(u!i); and from Eqs (3.7, 3.10 b) 



£bt(^i) = K - U - S B -$(ui)] 1 

= K - U}- 1 + K - C/r 2 S BT M + 0(11^) 

(3.13) 



where oj[ = u-y + irjSga(uJi — U/2). From Eq. (3.12), 
the leading large- U contribution to Q~^(lh) thus arises 



from the second term in Eq. (3.13), yielding Q B ^(u>) — 
U- 2 S^(uj); hence from Eqs ( |3.11|j3^ c): 



i \ a— i i ^ f 2**^At( w + w p) 



AF 
P 

(3.14) 



We focus now on the lower Hubbard band (LHB) 
in strong coupling, viz u> w <C U — > oo; the up- 
per Hubbard band follows trivially by symmetry. Since 
£ b t(w) = — S A |(?7 — u), it follows that for uj in the 
LHB £ B f(w) is pure real and 0(l/U); it can thus be 
neglected. The G^(ui » 0) are likewise pure real, with 
G+ T (w) - 0(l/U) and G+ T (w) - 0(l/U 3 ) (as may be 
shown using Eq. ( [3.10 ) together with the analogue of Eq. 
Q for together with G+(lo) = ±[G++G+], 

they too may be neglected. And from Eq. ( 3.10| b), 



G bT (w) = Q^(u>) = U S b ^(uj) which can also be ne- 



glected asymptotically. Hence in total, G A f(w) = G A Ju>) 
and G(oj) = i[G AT (w) + G bT M1 = \G^(uo). For the 
LHB in strong coupling, Eq. ( 3.10| a) thus reduces to 



G AT H = [u- ^^(uj + ujp)}- 1 
G At (w) = [w-^G AT (w)]" x 



AF 
P . 



(3.15a) 
(3.15b) 



These are the equations for the corresponding t-J z 
model on the Bethe lattice, for an AF and random spin 



7 



background respectively (see eg [|lT],^TJ ) ; the t-J z model 
itself is naturally equivalent for d°° to the t-J model since 
the spin excitations are purely Ising-like. We add in pass- 
ing that a much more detailed asymptotic analysis, pick- 
ing up constant terms Q (l/U) , leads to the 'bare' uj in 
the denominators of Eqs ( |3.15 ) being replaced by uj + uj„ 



and oj + ^ojp respectively for the AF and P phases. These 
shifts, neglected in the brief analysis above, reflect simply 
the presence of the trivial charge terms in the t-J model 
(see Eq. (^1])); they are irrelevant to our subsequent 
discussion. 

Since the t-J limit emerges correctly in strong cou- 
pling, the present theory is thus asymptotically exact. 
Consider for example the P phase, noting that for the 
U = non- interacting limit the Green function Go(u>) = 
ReGo(w) — iirsgn(oj)Do(oj) is given by 



G (oj) = [uj - i^Go(w)]- 1 : ^ = 
whence the non-interacting spectrum 



(3.16a) 



D (oj) = —[2 - (w/i*) 2 ]^ : M < V2U (3.16b) 
iris. 



is a semi- ellipse with full width 2 V2t*. From Eq. ( |3.15| b) 
this is also precisely the spectral density for G A j {uj) in the 
lower Hubbard band. And since G(oj ~ 0) = iG A j(w) as 
above, the total lower Hubbard band spectrum in strong 
coupling is L>l(w) = ^D (uj), see also §2. lb; (the nor- 
malization factor of \ naturally reflects the fact that 
the remaining half of the single-particle spectrum oc- 
curs in the upper Hubbard band centred on uj = U, viz 
Du(uj) = \D (U - uj)). Note further that the Feen- 
berg ('medium') and interaction self-energies contr ibute 
equally to the it 2 G A |(w) denominator in Eq. (3.15b) for 
the P phase. Physically, this reflects the fact discussed 
in §2. lb that, upon adding a er-spin electron to a site, 
it is equally probable for either the added er-spin or the 
— er-spin electron already present to hop off the site. At 
UHF/HIII level, in contrast, only the former canby con- 
struct occur (S QCT = 0): the analogue of Eq. ( 3.15| b) is 



then G^(uj) = [uj—jtlG^(uj)] 1 , producing an incorrect 
strong coupling bandwidth of 2t Sf as argued physically in 
§2.1b. 

The AF case itself is discussed further in the follow- 
ing section since, in contrast to the P phase, the ap- 
proach to strong coupling is subtle and physically reveal- 
ing. Here we simply add that (i) in contrast to the P 



phase, the \tlG A -\{uj + uj p ) denominator in Eq. ( [j.lii a) 



for the AF stems solely from the interaction self-energy 
£ A f(u;). Thus at UHF/HIII level atomic limit behaviour 
arises (incorrectly), viz G a -^(uj) — 1/uj, as argued phys- 
ically in §2. la. (ii) Alth ough obtained explicitly for the 
Bethe lattice, Eq. ( |3.15| a) holds equally for the hypercu- 
bic lattice in strong coupling. This is because retraceable 
paths, which by construct are the only self-energy paths 
for a Bethe lattice, are for the d°° hypercube also the 
only paths which restore the Neel spin configuration; see 
also [291. 



3.3 Results 



At finite U the basic self-consistency Eqs, ( |3.2p and 
(|3.8| a), are solved numerically. We consider first the AF 
phase. 



A. Antiferromagnet 



For U/t* — 10, Fig. 8 shows the resultant lower 
Hubbard band, Dl(uj) = 7r _1 ImG(w); from particle-hole 
symmetry the upper band follows by reflection about the 
Fermi level, Du(uj) = Dl(U — uj). 

For the same w p -value (Fig. 3a, inset), Fig. 8 shows 
also the corresponding t-J z limit spectrum from Eq. 
( |3. 15 a) . As is well known |llj the t-J z spectrum is dis- 
crete (and to illustrate relative intensities is thus shown 
with height proportional to integrated weight). Physi- 
cally, this reflects the fact that the hole is pinned by the 
string of spin-flips its motion creates, leading therefore to 
spatially localized single-particle excitations and hence a 
discrete spectrum; mathematically, it is reflected in con- 
vergence of the continued fraction implicit by iteration of 
Eq. (|3~T5|a). 

Although the U '/t* = 10 spectrum evidently bears a 
close resemblance to its t-J z counterpart, it is by con- 
trast continuous. This persists for any finite U: with 
increasing interaction strength the individual sub-bands 
in Dl(uj) centre ever closely on their t-J z counterparts, 
and their integrated spectral weights tend to those of the 
t-J z limit; but they retain a finite width, reflecting de- 
localization of the hole. The peculiarities of U — » oo 
are further evident in the t-J z model itself, Eq. (3.15i). 



For any uj p > th e t-J z spectrum is discrete, while for 



(as in Eq. (3.15b) for the P phase) the spectrum 



is continuous: the point uj p = thus corresponds to a 
transition from localized to extended single-particle exci- 
tations, and since uj p — > tl/U as U — » oo it is clear that 
U = oo is a singular point. 

While the physical mechanism leading to delocaliza- 
tion of the hole at any finite U is not of course inherent 
in the t-J z model Eq. ( 3.15| a) itself, it is readily inferred. 
Consider the Neel spin configuration and imagine remov- 
ing, say, an "["-spin electron from an A-type site, i. The 
nearest neighbours (NN) to any |-spin site all all |-spins. 
Hence to leading order in U — the t-J z limit — the hole 
initially moves via a NN J,-spin electron hopping onto 
site i, creating thereon a spin-flip (with an associated 
exchange energy penalty) ; and the subsequent motion of 
the hole via such a correlated sequence of alternating NN 
spin hops, in leaving behind a string of upturned spins, 
would by itself render the hole spatially confined. 

At large but finite U there is however a small but non- 
vanishing probability amplitude, of order tl/U, for an j- 
spin electron on a second NN site, also A-type, to hop to 
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site i via an intervening j.-spin site: the hole thus moves 
two lattice spacings, to the second NN A-type site. Un- 
like the " t-J z processes" above, this does not entail a 
spin-flip with concomitant exchange penalty: the hole 
moves freely. 

This mechanism evidently leads to hole derealization 
and, in tandem with the t-J z processes, produces the 
strong coupling spectrum. Its formal origins reside in 
the passage from Eq. (B.ICa) to Eq. ( |3.15| a) for the 
AF lower Hubbard band in strong coupling, where the 
Feenberg part of the self-energy S a ^(uj) = ^G B f(w) 
was neglected. As seen readily from the asymptotics 
of §3.2, the leading corrections to Im5 A |(w » 0) are 
lmS A i(uj) = (£*/2J7) 2 ImG A f (w). It is these that embody 
the derealization described above, and lead to spectral 
broadening (contributions to ReS^^u;) are 0(1/U) and 
lead simply to residual energy shifts). Further, note that 
since the energetic width of the spectral broadening is 
naturally the smallest energy scale in strong coupling, 
the principal effect on the 'bare' t-J z spectrum is a small 
resonant broadening of the individual t-J z lines. This 
is seen in Fig. 8, and becomes clearer still with further 
increasing U. 

To our knowledge, the above mechanism is the only one 
which can lead to hole derealization for the d°° AF in 
strong coupling; and for the reasons already given in §3.2 
applies to the hypercubic as well as the Bcthc lattice. In 
finite-d it is for example well known that Trugman paths 
p2[ lead to hole derealization for the hypercubic lat- 
tice, but such processes are 0(d~ A ) and do not therefore 
contribute in d°° 

As U is decreased, the spectra continue to exhibit es- 
sentially strong coupling behaviour down to modest in- 
teraction strengths of ~ 2 — 3, and can thus be 
understood quantitatively starting from the t-J z limit. 
This is shown in Ref. (see eg Fig. 3(b) therein). 

With further decreasing U however, the spectra evolve 
continuously to a weak coupling form that shows no trace 
of remnant t- J z -like behaviour. The spectral gap closes 
only in the non-interacting limit whence, correctly, the 
system is an AF insulator for all U > 0. The full spec- 
trum D(uj) = Dl + Djj is shown in Fig. 9 for U jt % = 1, 
together with the corresponding UHF spectrum to which 
(as one expects) it is qualitatively closer, although the 
single-particle gap A g is reduced to 0.42 of the UHF gap 
A = *7| Mo |. 

Two further renormalizations have been performed to 
check the veracity of the above results. First, note that 
although the Green f unctions have been obtained self- 
consistently via Eq.s J3.2| , 3.8 ), the single-particle prop- 
agators occurring i n th e RPA n~+ that enters the self- 
energy kernel Eq. (3^ a), are themselves bare UHF prop- 
agators; see Fig. 2. To ensure the theory is robust, we 
have thus additionally renormalized the single-particle 
lines entering n~+ in terms of both the (self-consistent) 
full Green functions G aa and the host Green functions 
Q a a- The results in either case differ only quantitatively, 
and at low U, from those just described; see also below. 



The second renormalization concerns the local moment 
which, in the calculations above, has been set to its 
UHF value \fMy\- In weak coupling, van Dongcn ^] has 
examined perturbatively the Neel temperature and the 
moment magnitude \fi\ (the order parameter) for the d°° 
hypercubic lattice, and has shown that even for U — > 0+ 
these are reduced by a factor q of order unity (q ~ 0.28 
p3| ) below their corresponding UHF values. The present 
theory is not of course perturbative (eg the emergence of 
the AF spin-flip scale is intrinsically non-perturbative) , 
but it is certainly closer in spirit to van Dongen to renor- 
malize the moment beyond UHF level. This is quantita- 
tively important at low U , and is achieved by requiring 
that be determined fully self-consistently via (cf Eq. 

(H) 



I Ml 



dQ [-Dat(^) - -°Ai(w)] 



(3.17) 



where D Aa is the full (as opposed to UHF) spectral den- 
sity. 

For illustration Fig. 9 shows the Bethe lattice spec- 
trum at U/t* = 1, obtained with both and n~+ 
renormalized (the latter in terms of the full Green func- 
tions). The gap A g is further diminished, the ratio 
g = Ag/A being ~ 0.15; while the local moment 
is likewise reduced below its UHF counterpart, such that 
rn = |/u|/|/io| ~ 0.39. It is not unfortunately feasible 
to obtain numerically accurate estimates of g and m as 
U — ► (since and A rapidly become exponentially 
small). But for U/t* — 1 the UHF moment itself is 
accurately represented by its asymptotic U —* limit, 
|/io| = 8-\/2exp[— nt*l \f2XJ — 1], so the above result for m 
may be reasonably close to its limiting value. 



B. Paramagnet 

To obtain correctly the strong coupling limit for cither 
phase is, as has been shown, fairly subtle. But in con- 
trast to the AF, the approach to strong coupling for the 
paramagnetic phase is not. Fig. 10 shows the full spec- 
trum D{Cj) = — 7r~ 1 sgn(cD)ImG(cD) (u) = oj — U/2) for the 
P phase at U/t* = 8, 6 and 4, comp ared to the strong 
coupling t-J z limit from Eq. (3.15b). For lift* = 8, 
the strong coupling limit has in practical terms been 
reached: the Hubbard bands are essentially symmetri- 
cally centred on w = ±Z7/2 respectively, with widths 
W 



Woo = 2-\/2£* and a band gap of A g ~ 

* ; even for U /£* = 6 the departure from the 
asymptotic spectrum is relatively minor. With further 
decreasing U , however, the individual bands become in- 
creasingly asymmetric; and the gap tends to zero more 
rapidly than , signalling the collapse of the insulating 
phase. This we now discuss, adding that throughout the 
insulating regime the local moments are well developed 
Qui > 0.95), as in Mott's conception of a Mott insulator 
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4. DESTRUCTION OF THE MOTT INSULATOR 

Fig. 11 shows the resultant band gap, A g {U), for the 
paramagnetic insulator as a function of U/t*. A g {U) is 
found to vanish continuously at a critical U c = 3.41t*. 
Detailed numerical analysis shows the corresponding ex- 
ponent to be unity, 



The divergence of \B\ controls additionally the lo w fr c- 
quency behaviour of ReG((D) = X{Cj). From Eq. (3.5c), 
X{Cj) = — X{— Cj), whence its leading low-w behaviour is 



X{Cj) = 7icj 







(4.5) 



(with 7i = -|7i|). From Eqs (3.1-3) and (3.6b), G(<D) 
may be written generally as 



A g (U)~(U-U c ) v 



v = 1 



(4.1) 



G(w) = i{[<D-4t2G(c5)-E AT (a;)]- 1 



(4.6) 



and we note that the width of the critical regime is 
quite narrow: the behaviour Eq. (4.1) is seen clearly for 
(U — U c ) < 0.05i*, corresponding to gaps A g (U) < O.li*. 

The continuous closure of the gap is intimately con- 
nected to the divergence of low-frequency dynamical 
characteristics of the system. Consider first the self- 
energy E A |((D). At frequencies Cj € [— £)+,£)+] inside the 
spectral gap (A a = 2cD + ), E A |(lj) = E A j(w) is pure real 
with a leading low-ci) expansion 



+ p-ii2G(£) + E AT (-£)]- 1 } . (4.7) 



Using Eqs (t 
|7i| to be related to \B\; the result is 



and (4.5) on either side of (|4.6|) enables 



7i 



l + l^l 

H 2 



A 2 



(4.8) 



E^ T (w) -A = BQ 



We find that A 2 > ^t 2 for all U > U c , whence the diver- 
gence of |B| as U — ► U c controls that of |7i|, 







(4.2) 



I7i I 







(4.9) 



Here, A = E aT (w = 0) (= -\U\fi\ + E aT (0), sec Eqs 
(3.1,8)), and is finite for all U > U c (see also below). We 
wish to find the behaviour of B — — \B\ as U — > U c . 

This is obtained by a scaling analysis. Defining y = 
Cj/Cj+, it is found that as the gap closes (<D+ — > 0), 
E A f {Cj) — A obeys the scaling form 



This is further confirmed by a scaling analysis of X {Cj) 
itself. In direct analogy to that for E A -|-(£>) above, X{Cj) 
is found to satisfy the scaling form 



X{Cj) 



w% x{y) 



(4.10) 



UJ + 



f(y) 



a 



(4.3) 



with x{y) = —x{—y), from which Eq. (4.9) in particular 
follows. For \y\ > 1 the corresponding spectrum D{Cj) 
likewise shows the same scaling form as expected, with 



with exponent a — ^; i.e 
close to U c , with correspondingly different gaps A g {U) 



for different values of U D{Cj) ~ Cj 1 / 2 [y-l} 1/2 = [Cj-Cj+Y' 2 for y - 1 close to the 



2Cj + {U), the w-dependent functions [E^(o)) — A]/ 



-1/2 



UJ 



- 



lower edge of the upper Hubbard band; this, combined 
with the spectral representation of |7i|, 



plotted in terms of y — Cj/Cj + , collapse to a 'universal' 
function f{y). Four points should be noted about the 
scaling behaviour, (i) Good scaling is found in practice 
for gaps A g < O.lt*, consistent with the critical regime 
found above for closure of the gap. (ii) The scaling is 
not confined to frequencies 5 « 1 well inside the spec- 
tral gap, but encompasses the region of non-zero spectral 
density {\y\ > 1), certainly up to \y\ ~ 2. (Similar scal- 
ing with a = \ naturally occurs for ImE A f(<D), as follows 
from Kramers-Kronig; see also below), (iii) In numerical 
terms the scaling analysis is sufficiently accurate to dis- 
tinguish readily between an exponent of a = „ and, e.g., 
a = |. (iv) The scaling function f{y) is a finite, well 
behaved function of y = Cj/ Cj + , with f{y) ~ y for y — > 
as is evident from Eq. ( [4.2|) . 

From Eqs (4.2) and (|4.2|) it follows immediately that 
151 



l7il=2 



dio 



D{Cj) 



(4.11) 



leads again to Eq. (4.9). 

It is instructive to compare the above results with those 
obtained from both the simple HIII approximation dis- 
cussed in §2. IB, and with the resonance broadening con- 
tributions |l3| additionally included, which we refer to 
as HIII'. For HIII, A = ~\U and_LB| = —the ap- 
proximation is purely static. Eq. ( fLq ) becomes a cubic 
for G{Cj), leading as is well kno wn to A g ~ {U — U c ) 3 ^ 2 
|||. As is clear from Eq. (UJ) with \B\ = 0, the tran- 



sition occurs when A 2 {U C ) = \t 2 , i.e. U c = \/2t*, and 
| 7 i| - (U-Uc)- 1 



-2/3 



--1/2 



The HIII' ap_ 
(ft 



Droximation 

can also be shown to be of the form Eq. ( |4.6| ), but with 
a w-dependent E A |(a)) given by E A -|-(a;) = — ^ + t 2 G{Cj) 
at low frequencies (which is sufficient to analyze the crit- 
ical behaviour); so that A — ~ ¥r and \B\ — t 2 \"f\\. Since 



\B\ 







or \B\ - {U - U c y 1/2 from Eq. (|t]) (which wc have 
confirmed by direct calculation of B = {dT, A i{Cj)/dCj) ). 



(4.4) E A |((I)) is a simple linear function of G(<D), Eq. ( |4.6| ) 
again becomes a cubic for G{Cj); an d, as for HIII, the 
gap exponent v = § {L3|. Eq. §J) with \B\ = ^| 7l | 
yields |7i| = {A 2 - \tl)/{A 2 - It*). The transition thus 
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i.e. U c -- 

{U-U c 



as is well 

-1 A — 2/3 



occurs when A 2 (U C ) = §ij, 
known |l5|]; and, again, |7i| ' 

Both HIII and HIII' are thus in the same universal- 
ity class, reflected more generally i n the fact that in ei- 
ther case scaling of the form Eq. (1.1C) can be shown 
to hold, but with an exponent of a = |. Gros has 
recently extended Hubbard's hierarchical equation of mo- 
tion decoupling scheme to higher order. The critical ex- 
ponents are unchanged from those of HIII/HIII'; and the 
value of U c itself is barely changed from its HIII' value of 
U c /t* ~ 2.45. From the above discussion it is apparent 
that the present theory belongs to a different universality 
class from that of HIII or its extensions. 

In direct analogy to the AF phase discussed in §3.3, 
we have tested the robustness of our results by further 
sclf-consistently renormalizing single-particle lines in the 
polarization propagator H~+(uj) = °IL~ ^'/(1 — U °n~+) 
that enters the self-energy kernel, Eq. ( |3.8| a). To illus- 
trate what this involves, consider renormalizing °n~+ 
(and hence n^+) in terms of the self-consistent host 
Green functions Q a a- T he res ultant °llj+(u;) is then 
given generally by Eq. ( 2.12J ), with the bare (UHF) 
Green functions G° Arr now replaced by Q Aa - For u) = in 
particular, Eq. (2.14) likewise holds, but with the bare 
Dj^ a (aj) replaced by the renormalized spectral densities 
Pa^w) = -7i" _1 sgn(w)Imt/ A(T (d;); ie 



u u n-+( w = o) 



(4.12) 



As discussed in §2,3 the key feature of the paramag- 
netic insulator is the zero- frequency spin-fli p sca le. To 
preserve this, the local moment \fx\ in Eq. ( 4.12| ) is it- 
self renormalized to ensure that at each step of the self- 
consistent iteration scheme U°U^ a (uj = 0) = 1 (and we 
note that throughout the entire insulating regime, the re- 
sultant moment is also self-consistent in the sense of 
Eq. (3.17) to < 1% accuracy). The results of this further 



renormalization are found to differ negligibly from those 
we have reported above. 

Finally, to demonstrate the importance of the ut s = 
spin-flip scale, we have eliminated it: bo th b y (a) ne- 
glecting its contribution to £ A j(iI>) in Eq. (|3.8| b), retain- 
ing only E| t T oncr (o)); and (b) replacing n~+ by °n~+ in 
the self-energy kernel Eq. (B.Sa). Results obtained from 
(a) and (b) are very similar, but differ qualitatively from 
those reported above. In particular, although the self- 
energy remains ^-dependent, the resultant critical be- 
haviour is found to be that of HIII/HIII' — the gap closes 
continuously, but with an exponent v = |. This points 
clearly to the necessity of including the lo s = spin-flip 
scale throughout the entire insulating phase: not only in 
achieving the correct strong coupling limit (as in §3.2), 
but also in describing the destruction of the insulating 
state. 



5. DISCUSSION 



We now discuss the present work, particularly in rela- 
tion to the iterated perturbation theory (IPT) approach 
use and application of which has been 
extensive | . Although our theory of the Mott-Hubbard 
insulating phases, with its explicit emphasis on local mo- 
ments, is conceptually and technically distinct from IPT, 
some general points of marked contrast are evident. 

For the antifcrromagnetic phase we have emphasized 
the importance of the lu p spin-flip scale, inclusion of 
which is necessary to obtain even qualitatively reason- 
able results throughout essentially the entire range of in- 
teraction strengths, and in particular to recover exact 
strong coupling asymptotics. However IPT does not ap- 
pear to capture the AF spin-flip scale, presumably be- 
cause it omits repeated particle-hole interactions of the 
sort shown in Fig. 2 (which, as in §2.2A, are required 
to pick up the spin-flip). This is seen, for example, from 
the known inability of IPT to describe correctly the U - 
dependence of the Neel temperature ||, particularly in 
the 'Heisenberg' regime. 

For the paramagnetic insulator, the results of §4 also 
disagree qualitatively with those obtained from IPT; see 
in particular Q and the review ||. Within IPT the 
paramagnetic insulating solution is found to break down 
dis ■continuously (at a critical U c i = 3.67t*, where the 
IPT gap A g (U cl ) ~ 0.3**), and | 7 i| (Eq. (4.10)) remains 
finite at the transition. 

The same authors fll2| have recently examined the insu- 
lator via exact diagonalization (ED) on clusters of n s = 3, 
5 and 7 sites, extrapolated to n s — * oo assuming l/n s 
scaling behaviour. The resultant data suggest a continu- 
ous closure of the gap at a U c %/t* — 3.04 ± 0.35 and are 
consistent with A g (U) ~ (U — U c ); see also ||. Further, 
and independently of the gap analysis, the behaviour of 
|7i| has also been examined by ED noting (see Eq. 
(4.10)) that a divergence in |7i| implies a continuous clo- 
sure of the gap: 1 / 1 Tx | is found to show good scaling 
behaviour, and to scale to zero when n s is extrapolated 
to oo. 

The present theory evidently agrees with the inferences 
drawn from ED. These concur with our predictions (§4) 
that the gap closes continuously and with an exponent 
v = 1, that |7i| diverges, and (less importantly) the value 
of U c itself; note moreover that the ED gap [0] is in 
rather good agreement with the present work over a wide 
U range. As described in §4, inclusion of the ui s = spin 
flip scale is central in describing the destruction of the 
Mott insulator. That IPT appears unreliable close to U c 
H thus suggests an incomplete inclusion of the effects 
of this spin scale — which cannot be entirely absent since 
IPT does give the correct strong coupling spectrum Hj - 
although in physical terms the origin of the spin-flip scale 
within IPT is not transparent. 

To conclude, we have developed in this paper a the- 
ory for the T — Mott-Hubbard insulating phases of the 
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d°° Hubbard model, encompassing both the antiferro- 
magnetic and paramagnetic insulators. The microscopic 
perspective it affords hinges on the importance of low- 
energy scales for insulating spin-flip excitations. Their 
existence is physically natural within the explicit local 
moment picture intrinsic to the theory, and inclusion of 
them is required not only to obtain the strong coupling 
limits of the single-particle spectra — which are captured 
exactly — but more generally to describe the entire insu- 
lating regimes, including for the paramagnetic phase in 
particular the destruction of the Mott insulator. 

Let us also note what we have not considered: the 
metallic state of the paramagnetic phase. But a glimpse 
of what is required to describe the metal within the 
present framework is evident from Fig. 12. For U/t* — 
3.5, close to the critical U c of §4, this shows the spectral 
density of transverse spin excitations Imn+~(cj) (here 
obtained, as described in §4, with °n+~ rcnormalized in 
terms of the Q Aa ). The w = spin- flip pole character- 
istic of the paramagnetic insulator is evident, and per- 
sists down to U c . Clearly, however, the spectral edges of 
the Stoner-like bands are themselves approaching uj = 0. 
This they do at U — U c , and for U < U c in the metallic 
phase the insulating spin-flip pole at u> = is replaced 
by a resonance at a small non-zero frequency u> = luk, 
indicative of the Kondo-like physics known to dominate 
the correlated metal ||||] • Extension of the present ap- 
proach to describe the metal, encompassing the Kondo 
spin-scale in such a manner that the correlated state is 
correctly a Fermi liquid, will be described in a subsequent 
paper. 
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FIG. 1. UHF single-particle spectrum, D°(uj), vs 
uj = uj-U/2 (in units oft,) for d°° Bethe lattice. At U/U = 4, 
for AF phase (full line) and P phase (dashed). 



FIG. 2. Particle-hole ladder sum in transverse spin chan- 
nel, for RPA n^~. Bare (UHF) propagators are denoted by 
solid lines, on-site interactions by wiggles. For d — oo, all 
intermediate sites i\ . . . i n are equal to i. 

FIG. 3. ImIL+A H vs u/U at U/U = 4 for d°° BL. (a) For 
AF phase; inset shows U/U dependence of AF spin-flip pole 
tup, with dotted line denoting U — * oo asymptote = tl/U. 
(b) For P phase, where spin-flip pole lu s = for all U in 
insulating state. 

FIG. 4. QMC Neel temperature vs U/U (open circles) 
for d°° hypercubic lattice |^,^5|. The simple estimate 
Tjv — |w p ({7), argued to be valid for U/t*>3, is also shown 
(solid line for U/t, > 3). The strong coupling asymptote 
t n = t*/ 2U is indicated by the dotted line. 



FIG. 5. Diagram contributing to single-particle self-energy 
Sio-, with same notation as Fig. 2; for full discussion, see text. 

FIG. 6. (a) Undressed (or self-consistent host) Green func- 
tion Qii-cr, expressed in terms of bare (UHF) propagators and 
site-diagonal self-energy insertions Sj CT . Note the restrictions 
on intermediate site sums, (b) Basic approximation to Ei CT 
used in present work, from which full set of diagrams retained 
follows by iteration using Fig. 6(a). 

FIG. 7. A 'direct' diagram (a), and its exchange counter- 
part (b). With site j = i excluded from the implicit sum over 
j, the direct diagram is O(l) while the exchange diagram is 
0(l/d). 

FIG. 8. Lower Hubbard band spectrum Dl(u)) vs ui (in 
units of i„) for AF phase (Bethe lattice) at U/t* = 10; the 
Fermi level lies at U/2 — 5. The corresponding t-J z limit 
spectrum is also shown, as discussed in text. 

FIG. 9. Full spectrum D(u>) vs uj — uj — U/2 (in units of t*) 
for AF phase (Bethe lattice) at U/U = 1. Dotted line: UHF 
spectrum; full line: from present theory; dashed line: with 
and nit further renormalized as described in text. 



FIG. 10. Full spectra D(u) vsii = w- U/2 (in units of U) 
for P phase (Bethe lattice) at U/U = 8 (a), 6 (b) and 4 (c). 
Corresponding strong coupling spectra are shown as dashed 
lines. 

FIG. 11. Resultant spectral gap A g (U) vs U/U for P 
insulator (Bethe lattice). The strong coupling asymptote 
= U — 2v2t* is also shown (dashed line) . 

FIG. 12. Imn+;(w) vs iv (in units of U) at U/U = 3.5 
close to the boundary of the P insulating state, with renor- 
malization as described in text. 
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